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ABSTRACT 


The  efficiency  oi  multilayer  analysis  in  calculating  resistivity  profiles  from 
spreading  resistance  measurements  depends  on  the  rapid  numerical  evalua¬ 
tion  of  the  well-known  correction  factor  integral  first  introduced  by  Schu¬ 
mann  and  Gardner.  We  present  a  new  approximate  form  for  the  correction 
factor  which  allows  its  numerical  evaluation  with  only  22  integrand  values 
for  each  evaluation  of  the  integral.  When  this  technique  is  used  in  our  multi¬ 
layer  analysis  program,  we  unfold  spreading  resistance  profiles  at  the  real 
^  time  rate  of  less  than  3  sec/point  on  a  desktop  computer.  Its  results  match 

^  those  of  analytic  evaluations  of  special  two  layer  cases  or  with  more  elaborate 

numerical  evaluations  of  graded  structures  to  within  1%. 


Two-probe  spreading  resistance  measurement  has 
become  a  widely  used  technique  for  determining  im¬ 
purity  profiles  in  semiconductors  (1-4).  The  conduc¬ 
tivity  of  the  material  contacted  by  the  probes  is  given 
by 

R  =  C/(2<ur)  [1] 


where  R  is  the  measured  spreading  resistance,  a  is  the 
contact  radius,  and  C  is  a  correction  factor  which  is  a 
function  of  the  conductivity  profile  of  the  material.  The 
multilayer  analysis  technique  of  evaluating  this  cor¬ 
rection  factor  was  originally  proposed  by  Shumann 
and  Gardner  (1)  and  later  developed  in  more  easily 
calculated  forms  (2-4) ;  complete  derivations  are  given 
in  Ref.  (2,  3).  In  this  technique  the  correction  factor 
U 


I  is  the  total  current  fiowing  between  the  probes, 
where  Ji  and  Jo  are  Bessel  functions  of  the  first  kind, 
D  and  t  are,  respectively,  the  probe  separation  and 
layer  thickness  measured  in  probe  radii,  F  (sometimes 
denoted  by  1  4-  2e)  is  a  function  of  all  conductivities 
deeper  than  the  contacted  layer,  and  f  (x)  is  the  Hankel 
transform  of  the  normal  component  of  the  current  den¬ 
sity  I(r)  at  a  probe.  F  can  be  expressed  by  a  recursion 
relation 

„  Fiiri+i  4- <ritanh(tx) 

Fu  I  = -  [3] 

»i  4-  Fivi+i  tanh(tx) 


*  Kteetroehamleal  Society  Active  Member. 

Key  words:  reilsUvlty  proflles,  spreading  resistance  measure¬ 
ments. 


which  is  equivalent  to  a  form  given  by  Choo  (3).  I(x) 
is  given  by 
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I(z)  =  Jo(.zr/a)I(r)rdr 


Two  forms  of  l(x)  are  found  in  the  literature 
I  Bin  lx) 


Hx)  =• 


2nx 


[4] 


which  was  first  used  by  Shumann  and  Gardner  (1)  and 


Ux)  =■ 


nx 


[5] 


which  was  suggested  by  Severin  (5)  and  used  by  Choo 
(6) ;  the  latter  is  used  in  the  remainder  of  this  paper. 
A  complete  discussion  of  the  consequences  of  using 
either  form  is  given  in  Ref.  (7). 

The  calculation  of  a  conductivity  profile  from  a  set 
of  measured  spreading  resistances  proceeds  as  follows. 
Starting  at  the  deepest  point,  the  value  of  the  conduc¬ 
tivity  (and  the  resulting  correction  factor)  is  adjusted 
until  the  calculated  spreading  resistance  agrees  with 
the  measured  spreading  resistance  to  the  required  ac¬ 
curacy.  When  this  is  achieved,  the  process  is  repeated 
at  the  next  point,  continuing  until  the  entire  set  of 
measurements  has  been  exhausted.  The  trial  conduc¬ 
tivity  enters  the  calculation  of  the  correction  factor 
only  through  the  function  F.  Hence  during  a  numerical 
evaluation  of  C  only  the  function  F  need  be  recalcu¬ 
lated,  the  remainder  of  the  integrand  in  Eq.  [2]  hav¬ 
ing  been  calculated  once  and  stored  as  a  numerical 
array. 

Because  the  multilayer  analysis  scheme  may  require 
several  evaluations  of  the  correction  factor  per  point 
an  efficient  technique  for  its  evaluation  is  desirable.  In 
this  paper  we  introduce  an  approximation  to  the  cor¬ 
rection  factor  integral  which  can  be  evaluated  rapidly. 
We  then  develop  a  numerical  method  which  evaluates 
the  approximate  integral  by  sampling  the  integrand 
at  only  22  points.  We  then  compare  this  method  with 
an  analytic  evaluation  for  special  cases  and  with  other 
numerical  methods  for  a  full  multilayer  case. 

Approximation  of  the  Correction  Factor 

The  correction  factor  integral  in  Eq.  [2]  contains  a 
term  proportional  to  Jo(.Dx)/2  which  describes  the 
voltage  reduction  at  one  probe  due  to  the  other  probe. 
On  physical  grounds  its  effect  can  be  signiftcant  only 
when  the  conductivity  of  the  contacted  layer  is  very 
high  or  when  the  spacing  between  probes  is  small,  for 
which  case  Choo  (3)  showed  that  the  primary  effect  of 
the  Jo  term  is  to  assure  that  the  integrand  vanishes  at 
*  =  0,  even  when  F  is  singular  at  the  origin. 

In  this  section  we  show  that  it  is  possible  to  ap¬ 
proximate  the  correction  factor  integral  by  an  integral 
without  a  Jo  term  by  integrating  from  a  nonzero  lower 
limit 

C.  =  J’J’  dx  )*  F(ti)  [«] 


FA  =  l/(tx).  Fn  =  tx 


allowing  the  integral  in  Eq.  [2]  to  be  evaluated  ana¬ 
lytically.  For  case  (a)  the  function  Fa  bes  a  pole  at 
X  =  0;  the  lower  limit  L  is  chosen  so  as  to  minimize 
the  difference  between  C  and  C.  for  this  case.  Equa¬ 
tion  L2J  written  as 


C 


8  lim  r 

“  Jl(X) 

r  Ji{x) 

Jo(Di) 

1  dx 

X 

X 

2 

J  x‘-« 

a  special  case  of  the  Weber  Schafheitlin  integral  (8), 
yields 


C=— (ln(D) +1/4)  [8] 

xt 


The  corresponding  value  of  C.  is  obtained  by  evaluat¬ 
ing  Eq.  [6] 


<-*o 


Ji^Cx) 


dx 


p  Jl*(x) 

Jo 


dx 


] 


which  yields 

2 

Ca  =  —  [ln(2e-VL)  +  1/4  +  WB . . .  ]  [9] 

xt 


where  y  is  Euler’s  constant,  0.577 _ 

With  L  set  at 

2e-y/D  =  1.12292/D 

the  fractional  difference  between  C  and  Ca  in  this  case 
is  minimized;  its  magnitude  can  be  estimated  by  evalu¬ 
ating  the  term  of  order  L*,  namely 

c  —  c 

- ^ - 0.14/ ( (In  D  +  0.25)D*)  [10] 

Using  this  value  of  L,  both  C  and  Ca  were  evaluated 
for  cases  (b)  and  (c);  results  are  shown  in  Table  I. 
These  three  cases  represent  very  different  physical 
situations.  The  fact  that  the  errors  introduced  by  using 
Ca  rather  than  C  are  smaller  than  1%  when  D  is 
greater  than  10  suggests  that  the  approximation  is 
generally  useful:  in  most  probe  arrangements  D  is 
greater  than  30.  The  numerical  procedure  for  evalu¬ 
ating  Co  is  much  simpler  than  that  necessary  to  evalu¬ 
ate  C.  The  latter  involves  calculating  the  integrand  at 
enough  points  to  follow  the  rapid  oscillations  of  the 
function  Jo(Dx).  For  example,  when  D  is  equal  to  30, 
Jo(Dx)  has  more  than  30  zeros  at  values  less  than  the 
first  zero  of  Ji(x). 


The  Numerical  Integrator 

In  this  section  we  describe  a  numerical  integrator 
which  has  been  tailored  to  take  advantage  of  the  prop¬ 
erties  of  the  integrand  of  [6]  in  order  to  minimize  the 
required  evaluations  of  F. 

Our  objective  is  to  choose  a  set  of  numbers  Xi  and 
weights  tOi  such  that  Ca  can  be  well  represented  by  a 
numerical  approximation. 


where  L  is  given  by 

L  =  2e-v/D  [7] 

Here  y  is  Euler’s  constant,  0.577... .  We  demonstrate  this 
by  direct  calculation  of  C  and  Ca  in  three  cases  which 
are  amenable  to  analytic  methods,  namely  (a)  a  con- 
c>.‘*ng  layer  over  an  insulating  substrate,  (b)  a  con¬ 
duce.  '  layer  over  a  perfect  conductor,  and  (c)  a  uni¬ 
form  semi-infinite  slab.  The  appropriate  forms  of  F 
are,  respectively 

FA  =  coth(tx),  FB  =  tanh(ti),  Fc  =  1 

For  any  multilayer  case  F  must  be  within  the  en¬ 
velope  coth(tx)  and  tanh(tx).  The  test  cases  therefore 
represent  the  extreme  cases  and  one  intermediate  case. 

For  thin  conducting  layers  (i.e.,  in  the  limi";  t-»  0), 
the  first  two  forms  can  be  approximated  by 


Cn  =  2roiF(txi)  [11] 

The  appropriate  numerical  integration  technique  was 
selected  so  as  to  accommodate  the  most  singular  case, 
F  =  coth  (ti)  taking  care  that  other  cases  were 
treated  with  sufficient  accuracy. 

Based  on  properties  of  (Ji/x)»,  we  divide  the  domain 
cf  integration  into  four  regions 


(A) 

L^x^0.1 

(B) 

O.l^x^l 

(C) 

1  ^  X  ^  jl.l 

(D) 

ii.i  <  X 

where  ji,i  is  the  ith  zero  of  the  Bessel  function  Ji  and 
ii.o  =  0,  ji.t  =  3.83171.  In  cases  of  practical  significance 


Vol  IZt,  So.  5 


MULTILAYER  ANALYSIS 


1U9 


Tgbit  I.  Functional  formi  of  eorroction  facton  with  L  =  1.12292/0 


F  C  C.  II  -  C./C  I 


1/xt 

1 

xt 


'  (In  O  +  0.2S) 


,  \  s,  4D  V  SO* 

4t/r 


a  /  \ 

- 1  In  D  +  0.25  +  -  +  ...  ) 

rt  '•  a  ' 

w  '  8»  4  ' 


0.14 


t(4  -  L»)/t 


l>>(In  O  4  0J8) 
0i)72/D 
0M/O> 


I 


i 


I 
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L  is  restricted  to  the  domain  0.001  ^  L  ^  0.1,  In  each 
region  we  choose  a  numerical  integrator  which  will 
assure  adequate  convergence  of  the  integrals. 

In  region  (A)  we  replace  (Ji/x)*  by  1/4.  In  this 
region  the  most  difficult  integrand  that  can  arise  will 
be  proportional  to  coth  (to.),  which  approaches  l/(tx) 
for  small  x.  To  treat  the  pole  at  x  =  0,  we  choose  In  x 
to  be  the  integration  variable;  this  segment  of  C,  is 
then 

—  I  ,  — F(tx)d(lnx) 

„  •/*=L  4 

We  then  apply  a  nine  point  Newton>Cotes  integrator  to 
the  integrand.  The  nine  points,  equally  spaced  in  In  x, 
are  given  by 


The  approximation  is  exact  if  F  varies  at  most  linearly 
with  X.  We  took  five  subintervals  between  successive 
zeros  of  Ji  to  h,«,  x  =  19.616,  and  found  that  the  in¬ 
tegral  truncated  at  this  point  was  adequate  for  most 
integrands.  [Albers  (10)  truncates  at  x  =  20  and  Choo 
(6)  truncates  at  x  =  20.8.]  However,  in  cases  where 
F  is  small  in  region  (A)  (for  example  when  F  = 
tanh  (tx)),  an  error  of  as  much  as  4%  can  be  made 
due  to  this  truncation.  By  using  one  additional  point 
at  the  centroid  of  (/t/x)^  taken  between  19.616  tmd 
infinity,  this  error  can  ^  reduced  to  below  1%. 

This  segment  of  Cn  is  then 

8 

—  2«)DjF(ti)  [17] 

x 


Xi  =  L  (0.1/L)W«>,  0^i^8 
This  segment  of  Cn  is 

s 

—  ln(0.1/L)  2  to*jXiF(tXi)  [12] 

4n  0 

The  weights  u)9.i  for  a  nine  point  Newton-Cotes  inte¬ 
grator  are  given  in  Ref.  (9). 

In  region  (B)  the  slowly  varying  (Ji/x)*  is  explic¬ 
itly  calculated.  Again  the  most  difficult  integrand  is 
proportional  to  l/(tx),  and  we  make  the  same  trans¬ 
formation  as  in  region  ( A) .  To  achieve  the  same  den¬ 
sity  of  points  in  (B)  as  in  (A),  with  L  =  0.001,  we 
use  a  five-point  Newton-Cotes  integration  formula. 
This  segment  of  Cn  is  then 

4 

-?-ln(10)  X  /i2(Xi)F(xi)  wyi/x,  [13] 
n  0 

The  weights  toj.i  are  given  in  Ref,  (9)  and  xt  = 
0.1  (10) 

In  region  (C),  (Ji/x)*  is  falling  rapidly;  a  five-point 
Newton-Cotes  integrator  in  x  rather  than  In  x  provides 
adequate  coverage.  This  segment  of  Cn  is 

s 

—  (jl.t  —  1)  2)  (Jl(Xl)/Xl)*  F(X|)108,i  [14] 

n  0 

Only  four  of  the  five  points  need  be  calculated,  as 
•It(ji,i)  is  by  definition  zero. 

In  region  (D),  the  rapid  decrease  in  (Ji/x)*  guaran¬ 
tees  that  the  integral  beyond  the  first  zero,  Ji.i,  will 
make  only  a  small  contribution  to  the  correction  factor. 
Using  the  mean  value  theorem 

F(x)  du(x)  =  F({)  (tt(b)  -  u(a))  [15] 

•fs»a 

where  a^i^b  and  du  =  (/i/x)*dx,  we  approximate 
the  integral  between  successive  zeros  of  Ji  by  choosing 
{ to  be  the  centroid  of  (Ji/x)*  between  zeros 

f  Jj*(x)di/x 

S 

f  =  — -  [16] 

(Jj(x)/x)*dx 


where  the  weights  tcn.!  and  ft  are  shown  in  Table  IL 
As  the  points  at  the  region  boundaries  need  be  calcu¬ 
lated  oidy  once,  the  complete  integration  requires  only 
22  evaluations  of  F. 

Testing  Procedures 

The  integration  technique  has  been  tested  in  the 
range  0.001  ^  t  ^  1  and  10  £  D  ^  1000  for  two  layer 
cases  where  analytic  results  can  be  obtained  and  for 
representative  multilayer  distributions.  We  present 
the  results  of  the  two  layer  calculation  and  a  sample 
multilayer  case. 

Test  1.— In  the  first  test  we  compared  Cn  calculated 
using  the  22  point  integrator  with  C,  calculated  ana¬ 
lytically  in  the  section  Approximation  of  the  Correc¬ 
tion  Factor. 

F(x)  is  replaced  by  tx,  1,  and  l/(tx)  in  the  cases 
now  under  consideration. 

We  evaluated  Cn  at  12  values  of  D  between  D  =  10 
and  D  =  1000  and  list  the  maximum  deviations  in 
Table  III.  Maximum  errors  were  found  on  comparing  C 
with  C,  and  Ca  with  Cn.  In  all  cases  the  maximum 
error  occurs  at  D  =  10. 

In  all  cases  we  find  agreement  between  the  approxi¬ 
mation  and  the  numerical  integrator  to  be  better  than 
0.3%. 


Table  II.  Centroidi  and  weights  used  in  region  (D) 


* 

q  =  jl.l 

b  as  jltUl 

Centroid 

t> 

Weight 
100. 1 

1 

3.83171 

7.01558 

5.28853 

0.688548E-2 

2 

7.01538 

10.1735 

8.480*8 

0.163S98E2 

3 

10.1735 

13.3237 

11.6681 

0.627777E-S 

4 

14.3237 

16.4706 

14.829 

0.8051S9E.8 

5 

16.4706 

19.6199 

17.9881 

0.173318E4 

6 

19.6159 

00 

39.3 

O.SIUE-S 

Table  III.  Maximum  deviations 


I  C  -  C.  I  I  c.  -  C.  I 

F  C  C. 


tx 

0.008 

0.008 

1 

0.007 

0.0008 

1/tx 

0.0006 

0.0008 
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Test  2.— To  measure  the  performance  of  our  22  point 
integrator  for  cases  where  analytic  techniques  were 
unavailable  we  constructed  a  Simpson  integrator  of 
1597  points  constructed  as  shown  in  Table  IV. 

This  integrator  was  programmed  in  three  versions 
to  accommodate  the  various  tests  needed: 

SIM  which  numerically  integrates 

SIM  =  —  1  ( - - - I - Fdx 

a  •'0  \  X  2  /  X 

SIMSG  which  evaluates  the  Schumann  and  Gardner 
correction  factor 

SMSG  =  i.  r 

Jt  VO  \  *  2  /  * 

SIMT  which  corrects  SIM  for  truncation  error  as 
described  in  the  section  The  Numerical  Integrator. 

Truncating  the  integral  at  z  =  20  allows  us  to 
compare  our  results  with  those  of  Albers  (10),  who 
generated  spreading  resistance  data  (MD  in  Table  V) 
from  model  resistivity  profiles  by  using  the  integration 
scheme  of  D’Avanzo  et  al.  (4)  over  the  interval  0  ^  z 

Table  IV.  The  1597  point  integrator 


Number 

Xt  min.  X,  max  Weights  of  points 


0 

0.001 

Trapexoldal 

0.01 

0.01 

Simpson 

0.01 

0.1 

Simpson 

0.1 

1 

Simpson 

1 

20 

Simpson 

Table  V.  Summary  of  result*  calculating  resistance  os  a  function 
of  depth  in  a  layered  medium  using  the  various  integrators 
mentioned  in  text.  All  calculations  ore  bosed  on  the  conductivity 
profile  plotted  in  Fig.  t  and  excerpted  in  this  table  for  the  cose 
0  =  2Mm,0  =  25,  f  =  0.015. 

A.  Comparison  of  our  large  Simpson  Integrators  SIH  and  SIMSG 
with  the  IntegraUon  scheme  of  D’Avanzo 

A. 


Alber*a  model  data 


Layer 

No. 

Condue- 

Uvtty 

(mbo<m) 

MD 

fpreadlng 

resistance 

(0) 

rraottonal  deviation 

1.SIHSG/HD 

I-SIM/MD 

1 

7.075E-03 

1.744E  +  04 

1.121E~06 

•-U53E--01 

4 

1.940E--0a 

3.420E  +  0S 

1.144E-03 

-1.0S4£v.oi 

7 

7.050S-01 

3.S9SE  +  02 

1.310E-04 

•'4.69I£«02 

10 

9.77dE-fOO 

3.044E  +  03 

1.780E>04 

-2.769E-02 

13 

4.000E4'01 

3.101E  +  0Z 

1.740E>04 

''2.534£'-02 

16 

1.010E4  02 

3.S30E  02 

1.924E-04 

<-2.470£-02 

10 

i-asiE-foa 

S.430E  +  03 

1.186E>04 

-2.47SE-02 

23 

S.547E4-01 

1331E  +  03 

9.730E-05 

-2.487£~02 

26 

3.040E4>01 

4.780E  +  0S 

1.031E-04 

-2.957E-02 

26 

6.203E<fO0 

2.737E  +  04 

6.695E-05 

-2.982E-02 

31 

S.aOE-01 

1.B40E  +  0S 

6.066E>06 

-3.994E-02 

34 

1.177E-03 

3.0Z2E  +  0S 

>2.8aiE-00 

-8.439E-02 

37 

7.926E~03 

3.072E  +  0S 

>9.837E-06 

»a499E-02 

40 

7.906E-'(I8 

3.072E  +  0S 

-a.389E'06 

-8.SOIE-02 

B.  Comparison  of  the  reaulU  of  our  22  point  Integrator  and  the 
Oauae  Laguerre  Integrator  with  SUIT  and  SIH,  reapectlvety 


B. 


snrr-oRMs 

1-8IM/S1HT 

Fraettonal  deviation 

l-(22)/SIHT  1-GL/SIM 

1 

1.976B>04 

lJ<8E-02 

3JSSE-03 

-S.S39E-02 

4 

3.071E>06 

3J32E-0a 

3.229E-03 

-0.799E-02 

7 

4.103E-f08 

e.247E-03 

a.461E  04 

-4.040E-02 

10 

3.130E4  0a 

e.002E-04 

-1.2691-03 

-3.S91E-02 

13 

3.100E4>02 

I.TOOE-Oe 

-1.300B-03 

-3.903X-02 

16 

3.730E4-Q6 

e.ao4E-os 

-1.457E-03 

-3.73ZE-0t 

19 

6.77ra:4>02 

3.421E-0S 

-1.S42E-03 

-3.781B-09 

22 

1J64B4-0S 

2.023E-0S 

-2.320E-03 

-3.790E-02 

26 

4.67SE4^03 

2.241E-0S 

-3.43SE-03 

-S.874E-02 

26 

2.61914- 04 

1.S30E-03 

-7.900E-04 

-4.409E-02 

31 

1.9611  -f  00 

6.1S7E-0S 

4.3S0E-03 

-S.807E-02 

34 

IJ79B4^00 

7.S3SX-04 

3.00SE-03 

-7.09SE-02 

37 

tJ66l4^06 

9.9a4E-04 

3.570E-0S 

-7.0S5E-02 

41 

Maot-vos 

«.9«tE-04 

3.S73B-03 

-7.09SE-02 

20.  Table  V  and  Fig.  1  show  one  such  data  set  calcu¬ 
lated  to  represent  the  conductivity  profile  of  a  silicon 
wafer  containing  a  Gaussian  implant.  Simulated 
spreading  resistances  were  calculated  from  these  con¬ 
ductivities  for  the  case  a  =  2  ^ni,  D  =  25,  t  =  0.015 
using  SIMSG.  We  compared  these  to  the  set  of  resist¬ 
ances  suppUed  by  Albers.  The  largest  fractional  de¬ 
viation  found  between  the  two  data  sets  was  less  than 
2  parts  in  10'*.  Spreading  resistances  were  then  calcu¬ 
lated  for  the  same  data  using  the  22  point  integrator, 
SIM,  and  SIMT.  We  find  that  agreement  between  the 
22  point  integrator  and  SIMT  to  be  better  than  5  parts 
in  10^.  This  represents  the  total  deviation  between  C 
and  Cs  for  this  case.  A  comparison  between  SIM  and 
SIMT  results  indicate  that  the  maximtun  relative  de¬ 
viation  due  to  truncating  at  z  =  20  for  this  case  is 
about  2%. 

The  Gauss  Loguerre  Integrator 

To  demonstrate  the  value  of  the  lower  limit  cutoS 
approximation,  we  repeated  the  tests  performed  on  the 
22  point  integrator  using  another  sparse  point  numeri¬ 
cal  integrator  the  Gauss  Iiaguerre  integrator  (GL  in 
Table  V)  of  Choo  et  al.  (11).  This  integrator  retains 
the  term  in  Ja(Dz). 

In  Fig.  2  we  plot  the  fractional  deviation  between  C 
and  the  Gauss  Laguerre  integrated  correction  factors 
in  the  cases  of  the  section  Testing  Procedures,  Test  1. 

We  find  that  the  deviations  for  these  cases  tend  to 
cluster  in  a  band  of  ±  5%  with  occasional  excursions 
to  10%  or  20%.  This  is  due  principally  to  the  fact  that 
for  O  greater  than  10,  Jo(Dx)  has  more  than  60  zeros 
in  the  interval  of  0  ^  z  ^  20.  The  33  points  used  for 
the  GL  integrator  cannot  hope  to  cover  it  adequately. 

Discussion  and  Conclusions 

We  have  presented  a  new  approximate  form  for  the 
correction  factor  integral  of  Schumann  and  Gardner 
which  facilitates  the  rapid  and  accurate  evaluation  of 


Fig.  I.  Display!  data  supplied  to  us  by  Albers  (10) 

Fig.  lA.  Plot  of  condueHvily  vs.  depth  for  o  model  Goussien 
implant. 
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Fig.  1 1.  Flat  of  ipraading  vs.  depth  calculated  by  Albeis  using  the 
integration  scheme  of  D'Avanzo  (4).  The  set  shown  here  was  col- 
euloted  for  o  =  2  /m,  D  ^  2S,  t  z:  0.015. 


spreading  resistance  correction  factors.  This  has  been 
accontplished  by  eliminating  the  rapidly  varying 
Jo(Dx)  term  from  the  usual  correction  factor  integral, 
Eq.  [2]. 

Although  we  have  presented  the  approximate  in¬ 
tegral  of  Eq  [2]  only  for  I  (x)  satisfying 


/(*)  = 


IJt(x) 

It® 


the  approximation  remains  vaiid  with  the  lower  limit, 
L,  given  by  Eq.  [7]  for  all  current  density  profiles 
whose  Hankel  transform  satisfies  Eq.  [4b]  of  (5) 


/(*)  = 


lr(v+  l)Jv<.x) 

2ii(®/2)‘' 


All  of  the  commonly  used  forms  of  I(®)  are  special 
cases  of  the  above. 

The  integrator  of  the  section  The  Numerical  In¬ 
tegrator  was  incorporated  into  a  multilayer  analysis 
program  described  elsewhere  (12).  The  program  was 
written  in  BASIC  for  an  HP  system  9845.  The  time  to 


Fig.  2.  CorracHon  foctars,  GL,  were  colculoted  using  the  Gouts 
Laguerre  integrater  (ID  for  the  Shumann  and  Gardner  functian,  F, 
replaced  by  tx,  I,  1/tx  and  10  ^  0  1000.  Figure  2  shows  a  plat 

of  the  fractional  deviation  of  GL  from  the  analytically  calculated 
correction  factors,  C  of  Toble  I.  The  squores,  triangles,  and  circlet 
refer  to  cates  where  F  it  set  equol  to  tx,  1,  1/(tx),  respectively. 


calculate  one  point  was  generally  less  than  3  sec.  As 
the  time  to  measure  a  spreading  resistance  on  modern 
equipment  (e.g.,  an  ASR-100  probe)  is  about  6  sec  the 
complete  conductivity  profile  can  be  corrected  as  mea¬ 
surements  are  made. 

Manuscript  submitted  April  14,  1980;  revised  manu¬ 
script  received  Nov.  6, 1980. 

Any  discussion  of  this  paper  will  appear  in  a  Dis¬ 
cussion  Section  to  be  published  in  the  December  1981 
Journal.  All  discussions  for  the  December  1981  Dis¬ 
cussion  Section  should  be  submitted  by  Aug.  1,  1981. 

Publication  costs  of  this  article  were  assisted  by  the 
U.S.  Army  Electronics  Technology  and  Devices  Labora¬ 
tory  (ERADCOM). 
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